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Abstract. 

^ ' Introducing a surface layer of matter on the edge of a neutron star in slow rigid rotation, 

we analyze, from an intrinsic point of view, the junction conditions that must be satisfied 
^ ' between the interior and exterior solutions of the Einstein equations. In our model the core- 

{^jq' crust transition pressure arise as an essential parameter in the description of a configuration. 

As an application of this formalism, we describe giant glitches of the Vela pulsar as a result of 
variations in the transition pressure, finding that these small changes are compatible with the 
^ ^ expected temperature variations of the inner crust during glitch time 

m 

1. Introduction 

Many pulsars show sudden spin jumps, glitches, superimposed to tlie gradual spin down due to 
the continued loss of angular momentum suffered by the star. The study of the properties of 
the star during the glitch time is essential to understand the structure of the neutron star that 
models the pulsar [ll|2]. In Vela pulsar, giant glitches with relative period variation of the order 
; of 10-6 have been observed. Many mechanisms have been proposed to explain glitches. 

^ ' The models that deal with glitches are associated with the layer structure of neutron stars. 

^ Two regions of the star can be differentiated: the core and the crust. It is thought that the 

crust region has a solid crystalline structure similar to a metal [Bj. Some theories propose that 
the glitch is triggered by the rupture of the crust as a consequence of the tensions on the crust 
that try to adequate the ellipticity of the crust to the changing angular velocity of the star. 

Because of the large density gradient close to the transition between the core and the crust, 
most of the crustal matter resides in the shells in the inner part of the crust [4:J. Furthermore, it is 
noteworthy that the dynamical properties of the neutron star depend strongly on the transition 
pressure between the star's core and the crust as have been pointed up by [5j and [6j. Because 
most of the matter of the crust is found in the shells near the transition region between the core 
of the star and the inner crust, this is the region where the properties of the crust are relevant 
[3]. Hence, we will treat the crust as a surface layer that envelops the star's core. In the next 
section we will describe how our model of slowly rotating neutron star with a surface layer crust 
is constructed. 



2. Construction of models of slowly rotating neutron stars with surface crust. 

We assume that the neutron star is in permanent rigid rotation. The star rotates with constant 
angular velocity O around an axis, so the resulting space-time is axisymmetric and stationary. 
We assume that the rotation of the star is sufficiently slow and, hence, the Hartle-Thorne 
perturbative solution can be used [7, 8j. To second order in the angular velocity of the 
star f2, appropriate coordinates can be chosen so that the metric can be written as follows 
= -e2'^W[l+2/i(r,0)](it2+e2A(r)[i^2m(r,0)]dr2+r2[l+2A;(r,e)] [d9^ - sm'^ 9 {d^ - uj{r)dt)^]. 
We have to solve the Einstein equations in the interior of the star, where the matter is 
described by a perfect fluid tensor with equation of state p = p{p) and in the outer region, 
where space is empty. Also, we have to impose appropriate boundary conditions in the surface 
of the star. It is very desirable to introduce a new radial coordinate R adapted to the surface 
of the star [71[H]. 

We assume that the core of the star is surrounded by a thin surface layer (simplified crust 
model). Then we find that in the surface of the star R = a, the stress-energy tensor is non 
null and it can be written as T>^'^{R) = pc{0)u'^u'^{a), i.e. the surface layer can be described 
as a perfect fluid with its own angular velocity ilc- The surface energy density can be written 
up to second order as Pc(d) — E -\- 5£{6^. We must impose the usual junction conditions on the 
interior and exterior solutions, but taking into account the introduction of a surface layer of 
matter on the border of the star. We use the intrinsic formulation of these conditions [9] |10| : 
the first fundamental form must be continuous across the surface of the star, and the second 
fundamental form presents a discontinuity given by the stress-energy tensor in the surface of the 
star. In order to use these junction conditions in the slow rotation approximation, we expand 
them in terms of il, finding the following results for each one of the metric functions: 

Order zero: we obtain important conditions for the mass function and for the pressure of 
transition between core and crust 
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The equation ([T]) gives the total mass of the star to zero order [Mext) as the addition of three 
terms: the core mass {Mint = 4'/r Jq pR'^dR), the crust mass, and a negative bounding energy 
term. We interpret ([2D as an equation giving the surface density of energy e in terms of the total 
mass of the star, the interior mass, and the core- crust transition pressure Pint- This condition 
is used to obtain the radius of the star. 

First order: we obtain continuity of the rate of rotation of the inertial frames (x'(o)ext — 
oj[a)int = oj{a), and an expression for the total angular momentum: 
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where uJ^ = Oc — a;(a). It can be seen 0!8] that the first term of expression ([3]) has the same sign 
as O (in our case always positive). However, the second term depends on the sign of Oc- Hence, 
contra-rotating configuration are possible (the core and the crust contra-rotate), and even, we 
could find certain critical configuration in which the total angular momentum is null. 

Second order: we obtain A[r*(a, ^)] = where r* is defined as in [318], so both the mean 
radius and the eccentricity are continuous. Also, conditions for the second order perturbation 
of the mass function and surface density {5e) are obtained. 

Finally the last of the second order conditions fixes the value of the angular velocity of the 
crust Q.C = u;(a) it (Q — a;(a)). We obtain two possible configurations for the same star core 
and transition pressure: one with the crust co-rotating with the interior fluid and with identical 
velocity, and a special configuration with contra-rotating crust. 



3. Some properties of the model using realistic equations of state for the core 
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Figure 1. Mass vs radius for static/rotating 
configurations with a 10% crust (BPS) 
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Figure 2. Angular momentum vs transition 
pressure for fixed pc (BPS). 



To obtain a model of a neutron star with a surface layer crust in permanent slow rotation, for 
a given equation of state of the core, we have to fix three parameters. We choose the total mass 
of the star, the central density, and the cove- crust transition pressure. Then, we integrate the 
Einstein equations inside the star (using a realistic equation of state for the core) and we match 
this solution with the exterior solution found in |71 |8j using the junction conditions explained 
above. We have performed our simulations using an equation of state for high pressures from 
[11, K=240 MeV,^ = 0.78 and = 0.6] and the BPS equation of state for the lower pressures. 
Similar results have been obtained for the analytical fit of the Sly equation presented in [3]. 

We present some of our results in figures [1] and [2j In fig. [H we show multiple configurations 
for both static stars and rotating stars, with a crust of approximately 10% of the total mass and 
mass shedding angular velocity. On the static configurations without crust we write the value 
log lO(pc)- Note that the mass-radius behavior of the configurations with surface crust resembles 
those of the quark stars [llj . In fig. [2] we show the angular momentum for mass shedding angular 
velocity as a function of the core- crust transition pressure. In the counter rotating configuration 
we find that if the core- crust transition pressure is large enough, the angular momentum of the 
inner crust region could be equal or even higher than the angular momentum of the star core, 
resulting in a null or negative total angular momentum. 

4. Core- crust transition pressure evolution in giant glitches of the Vela pulsar 

On a typical pulsar, the rupture of the solid crust produces a release of energy 10^^-10^'^ ergs to 
the inner part of the crust. Several works study the thermal response of a neutron star after a 
glitch [121 |13j , obtaining that increases in the temperature of the transition region of the order 
of W^K. This small rise of the temperature will produce an small but sensitive increase in the 
transition pressure [13]. We propose that the glitch can be explained by changes in the equation 
of the state of the star in the transition region between the crust and the core. In our model, 
these changes are represented by variations of the transition pressure due to the increase of the 
temperature produced by the deposition of the energy, coming from the rupture of the crust. 

Using the observational data from the ATNF Pulsar Catalogue [15] we can study the giant 
glitches of Vela pulsar using our model. We assume that during the glitch neither the total mass 
nor the central density of the star are changed. Hence, in our description, the changes in the 
pulsar during this epoch are due to readjustment of the cove- crust transition pressure, i.e., a 
variation of the equation of state in transition region between the core and the crust. Assuming 
this, we can give a description of the effects of these glitches over the neutron star using our 
model of slowly rotating neutron star with superficial crust. We assume that the star has a total 
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Table 1. Date of each giant glitch of the Vela pulsar in modified Julian days and relative 
increases of angular velocity, transition pressure, eccentricity and quadrupolar moment 

mass of IAAMq, a central density of 1.279 • 10^^ gcm~^ , and a transition pressure before the 
glitch of 3.751 • W^^ dyn / cm'^ . We calculate the resulting configuration after the glitch varying 
the cove- crust transition pressure until we reach an equivalent configuration (same mass and 
central density) with the resulting angular velocity just after the glitch. As it can be seen in 
table dl our model predicts a relative transition pressure increase of the order of 10~^^ for the 
considered glitches. This small variation in the equation of state on the surface of the star causes 
a redistribution of other properties of the star. The eccentricity of the star increases order 10~^. 
Also, a movement of matter on the surface of the star, from the poles to the equator, produce 
an increase of the quadrupolar momentum (10~^) which may cause gravitational radiation. 

Using a estimation for the thermal pressure variations it can be seen that the corresponding 
increase of the inner crust temperature is order W^K |14] . which agree with the results obtained 
in the works previously commented. So the relative core- crust transition pressure changes we 
obtain in our model can be understood as variations in the equation of state of the transition 
region between the core and the crust of the neutron star, resulting from the deposition of energy 
of Iff^^-lO^^erg in the outer layers of the neutron star during the glitch time. 

In conclusion, we propose a new model for the giant glitches of the Vela pulsar: The tensions 
on the crust trying to adequate the ellipticity of the crust to the changed angular velocity 
produce the rupture of the crust. Then, there is a sudden energy deposition in the inner crust 
layers, which causes a rise of its temperature, resulting in an increase of the cove- crust transition 
pressure of the neutron star (a change in the equation of the state in the cove- crust transition 
region is produced). Hence, the properties of the neutron star (angular velocity, eccentricity, 
surface matter distribution and quadrupolar moment) change and a glitch is generated/observed. 
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